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Abstract 

We present spatially-resolved properties of molecular gas and dust in a gravitationally-lensed 
submillimeter galaxy H-ATLAS J090311.6-I-003906 (SDP.81) at z = 3.042 revealed by the 
Atacama Large Millimeter/submillimeter Array (ALMA). We identified 14 molecular clumps in 
the CO(5-4) line data. The surface density of molecular gas (Shs) and star-formation rate 
(Ssfr) of the clumps are more than three orders of magnitude higher than those found in local 
spiral galaxies. The clumps are placed in the ‘burst’ sequence in the Srs-TIsfr plane, suggest¬ 
ing that 2 ; ~ 3 molecular clumps follow the star-formation law derived for local starburst galaxies. 
With our gravitational lens model, the positions in the source plane are derived for the molecu¬ 
lar clumps, dust clumps, and stellar components identified in the Hubble Space Telescope image. 
The molecular and dust clumps are confined within a ~2 kpc region, while the spatial extent of 
the stellar components is as large as ~6 kpc and offset toward the west. The molecular clumps 
have a systematic velocity gradient in the north-south direction, which may indicate a rotating 
gas disk. One possible scenario is that the components of molecular gas, dust, and stars are 
distributed in a several-kpc scale rotating disk, and the stellar emission is heavily obscured by 
dust in the central star-forming region. Alternatively, SDR81 can be explained by a merging 
system, where dusty starbursts occur in the region where the two galaxies collide, surrounded 
by tidal features traced in the stellar components. 

Key words: gravitational lensing: strong — galaxies: formation — galaxies: individual (SDP.81) — galax¬ 
ies: ISM — galaxies: starburst — submillimeter: galaxies 
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1 Introduction 

Deep and wide-field submillimeter surveys have uncovered 
dust-obscured star-forming galaxies at high redshift (referred 
to as submillimeter galaxies; SMGs) (e.g., Blain et al. 2002, for 
a review). SMGs have intense star-forming activity with star- 
formation rates (SFRs) of a few 100-1000 Mq yr“^, which 
could be triggered by gas-rich galaxy mergers (e.g., Tacconi 
et al. 2006). SMGs and local ultra-luminous infrared galax¬ 
ies (ULIRGs) are known to be on the same relation on the 
molecular gas mass-far infrared (FIR) luminosity plane (Mh 2 - 
^fir) or molecular gas surface density-SFR surface density 
plane (Ehs-^^sfr), so-called ‘burst sequence’, high above that 
of ‘normal’ star-forming galaxies at local and z ~ 1-3 universe 
(Daddi et al. 2010; Genzel et al. 2010). Spatially resolved ob¬ 
servations of molecular gas in SMGs are essential to investigate 
the star-forming properties. However, only a handful of studies 
exists (e.g., Swinbank et al. 2010; Sharon et al. 2013; Rawle et 
al. 2014; Hodge et al. 2015) due to the limited spatial resolution 
and sensitivity of existing instruments. 

Gravitational lensing is a powerful tool to probe the physical 
properties of distant galaxies in detail because of its magnifi¬ 
cation effect. The gravitational lensing caused by a foreground 
massive galaxy or cluster of galaxies has been used to probe 
dusty star-forming galaxies at high-redshifts (e.g., Smail, Ivison 
& Blain 1997; Smail et al. 2002; Negrello et al. 2010; Vieira et 
al. 2013). Swinbank et al. (2010) discovered a z = 2.3 SMG 
(SMM J2135—0102) which is lensed by a z = 0.3 foreground 
galaxy cluster with the magnification factor of /x = 32.5, allow¬ 
ing them to examine the properties of the SMG at 100 pc scale. 

H-ATLAS J090311.6+003906 (SDP.81) is one of the bright¬ 
est sources discovered in the Herschel Astrophysical Terahertz 
Large Area Survey (HATLAS) (Bales et al. 2010). Multi¬ 
wavelength follow-up observations have shown that SDP.81 is 
an SMG at z = 3.042 ±0.001 that is gravitationally lensed by 
a foreground massive elliptical galaxy at z = 0.2999 ± 0.0002 
(Negrello et al. 2010; Negrello et al. 2014). Corrected for the 
gravitational lensing magnification factor of /x = 11 (Bussmann 
et al. 2013), SDP.81 has an intrinsic FIR luminosity of Lfir = 
5 X 10^^ Lq and a SFR of ~500 M© yr“\ Recently, the 
Atacama Large Millimeter/submillimeter Array (ALMA) ob¬ 
served SDP.81 at 1-2 mm in the Science Verification observa¬ 
tions with an angular resolution up to 0."03 (which corresponds 
to ~200 pc at z = 3.042), offering an excellent opportunity to 
study the spatially-resolved properties of an SMG. 

In this Letter, we present the properties of SDP.81 re¬ 
vealed by the high-angular resolution observations with ALMA. 
Throughout the paper, we adopt a cosmology with Hq —12 km 
s“^ Mpc“^, Qm = 0.26, and Qa = 0.74, and 1" corresponds to 
7.78 kpc at z = 3.042. 


2 Data 

2.1 ALMA Data 

We used the calibrated data set publicly available through the 
ALMA Science Portal. The details of observations and data 
reduction are described in ALMA Partnership et al. (2015a) and 
ALMA Partnership et al. (2015b). SDP.81 was observed as part 
of Science Verification of the ALMA Long Baseline Campaign 
in October, 2014, at Band 4 (151 GHz, or 2.0 mm). Band 6 
(236 GHz, or 1.3 mm), and Band 7 (290 GHz, or 1.0 mm) with 
the baseline lengths of ~5-15 km and with 31-36 antennas. 

The data were reduced with the Common Astronomy 
Software Applications (CASA) (McMullin et al. 2007) pack¬ 
age. The maps were processed with the CLEAN algorithm with 
the Briggs weighting (with robust parameter of 1.0). In this 
Letter, we used the Band 4 CO(5-4) line and the Band 7 
1.0 mm continuum images without tapering the uv-data at long 
baselines. The achieved synthesized beamsize (full-width at 
half maximum) is 07060 x 07054 and 07031 x 07023 for the 
CO(5-4) and 1.0 mm images, respectively. The images are 
shown in Figure 1. Two arc structures in east and west caused 
by gravitational lensing are clearly seen. 

2.2 HST Data 

We made use of the archival data of HSTAVide Field Camera 3 
(WFC3) taken by Negrello et al. (2014). We retrieved the cal¬ 
ibrated, fiat-fielded individual exposure frames of WFC3 with 
the F160W filter (i7-band) from the archive, and then coadded 
them by ourselves using the IRAF MultiDrizzle package. 
The final pixel scale was set to 07064 to improve the resolu¬ 
tion and quality of the image. We calibrated the coordinates 
of the image using the catalog data from the Data Release 12 
of the Sloan Digital Sky Survey (SDSS)^ We used 46 SDSS 
sources located in the image to derive the astrometric fit using 
IRAF. The rms error to the fit is 0!'085 in right ascension and 
07l00 in declination. We checked that the coordinates of the 
foreground lens galaxy in the SDSS catalog match the coordi¬ 
nates of the continuum emission detected in Band 4, 6, and 7 
(ALMA Partnership et al. 2015b) within 0f'015. The offset be¬ 
tween the center of the foreground lens and the ALMA contin¬ 
uum emission is smaller than the ALMA beamsize, showing the 
consistency of astrometry between the HST and ALMA images 
(Figure Ic). 

We used the software GALFIT (Peng et al. 2002) to deblend 
the background lensed source from the foreground lens galaxy. 
According to Negrello et al. (2014), two Sersic profiles were si¬ 
multaneously fitted to the foreground lens galaxy to cleanly sub¬ 
tract it. As a result, the brighter component of the foreground 
lens galaxy has the Sersic index of n = 4.42 (elliptical), while 
the fainter component has the index of n = 1.03 (exponential). 
The profiles are consistent with those obtained by Negrello et al. 
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Fig. 1. (a) ALMA CO(5-4) velocity-integrated intensity image. Molecular clumps identified in this study are presented. Synthesized beam is shown at the 
bottom left. Magnified figures (l."2 x l."2) of (a-1), (a-2), and (a-3) are presented, (b) ALMA 1.0 mm continuum image. Dust clumps identified in Tamura et al. 
(2015) are shown as plus marks. Synthesized beam is shown at the bottom left, (c) HST/WFC3 1.6 /xm. Plus mark represents the position of the foreground 
lens galaxy detected in the ALMA 2.0 mm continuum image (ALMA Partnership et al. 2015b). (d) H5T/WFC3 1.6 fim lens-subtracted image with contours of 
the ALMA 1.0 mm image (from +3cr with 2a step). Plus mark represents the position of the foreground lens galaxy detected in the ALMA 2.0 mm continuum 
image . Crosses represent the positions of stellar components which we used for deriving source-plane positions. 


(2014). The stellar arc structure seen in the lens-subtracted im¬ 
age shown in Figure Id is similar to the arcs seen in the ALMA 
CO line and continuum images, but they are clearly offset to¬ 
ward the west. 

3 Lens Model 

The analysis of gravitational lensing was conducted by Tamura 
et al. (2015) using the GLAFIC software package (Oguri 2010). 
The lens model was constructed by identifying the knots which 
were common in all three continuum images at 1.0, 1.3 and 
2.0 mm. The lensed images are explained by the best-fit mass 
model with a mass profile of an isothermal ellipsoid with a 400- 
pc core. The core radius of the isothermal ellipsoid agrees well 
with that of the stellar light in the HSTAVFC3 image, and the 
position of the predicted mass centroid is consistent with that of 
the emission detected in the ALMA continuum images. Tamura 
et al. (2015) found that the continuum emission has a power- 
law spectral index of —0.64 zb 0.1 from 1.0 mm to 20 cm, and 

^ http://skyserver.sdss.org/dr12/ 


it is likely that the non-thermal emission arises from an active 
galactic nucleus of the lens galaxy. This lens model is used to 
compute the positions on the source plane (Section 5.2). 

4 Molecular Clumps 

4.1 Clump Identification 

We identified molecular clumps by using the line cube (spec¬ 
tral resolution of 21 km s“^) and the velocity-integrated inten¬ 
sity map of CO(5-4). We first searched for the peak positions 
with >3.5cr in the velocity-integrated intensity map created by 
including pixels above L5cr in the line cube. We then used 
the CLUMPFIND package (Williams et al. 1994) to search for 
clumps in the line cube. We identified 14 molecular clumps 
along the arcs detected by both methods (Figure la). The po¬ 
sitions, magnification factors, and mean velocities are shown 
in Table 1. We caution here that the gravitational lensing can 
possibly produce two (or multiple) clumps in the image plane 
that are in fact the same clump in the source plane. However, 
because of the low SN, it may be difficult to identify all of the 
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multiples in the CO(5-4) image. We use the positions in image 
plane and source plane, and the mean velocities, as a guide to 
identify multiple images, and we suggest that clumps 1 and 13 
and clumps 4 and 5 are possibly counter images of the same 
source. 

If the clumps are spatially resolved with the ALMA beam, 
the transverse direction should be more extended than the radial 
direction. We found that some of the clumps do not appear to be 
extended in the transverse direction, and it is likely that they are 
unresolved. In order to measure the size of the clumps, higher 
angular resolution observations are needed. 

The magnification of the four clumps close to the critical 
curve (ID 4, 5, 7, and 8) is large (>100) and varies significantly 
within each clump. It is possible that they have properties with 
different physical scales compared to the other clumps, and we 
exclude these clumps in the discussion of star formation prop¬ 
erties. Except for the four clumps, the intrinsic spatial scale is 
roughly a few 100 pc in the source plane after correcting for the 
magnification factor and assuming that the clumps are spatially 
resolved. 

4.2 Surface densities of molecular gas and SFR 

We derived the surface density of the molecular gas and the SFR 
in order to investigate the intrinsic properties of SDR81. We as¬ 
sume that there is no significant difference in the distribution 
of CO and dust emission within the measured region for each 
clump, and thus the effect of differential magnification (e.g., 
Blain 1999; Serjeant 2012; Hezaveh et al. 2012) is negligible. 
Because the distribution of molecular clouds or star-forming re¬ 
gions smaller than the beamsize is unknown, we assume a filling 
factor of unity in this study. We note that the filling factor might 
be less than unity if they are not spatially resolved. 

We used regions above 2a in the CO(5-4) velocity- 
integrated map for conducting photometry and estimating the 
physical quantities of each clump. The photometry of dust con¬ 
tinuum was conducted on the 1.0 mm image in the same regions 
used for the CO photometry. The CO luminosity of each clump 
was calculated from Lqq = 3.25 x 10 ^+ 
z)~^ (Solomon & Vanden Bout 2005), where Sco^v is the 
velocity-integrated intensity in Jy km s“^, z^obs is the ob¬ 
served line frequency in GHz, and Dl is the luminosity dis¬ 
tance in Mpc. The error of Sco^v is calculated by using 
the equation of Mainline et al. (2004). We used a CO line 
ratio of CO(5-4)/CO(l-0) =0.3 (ALMA Partnership et al. 
2015b). The molecular gas mass is derived from Mh 2 = 
acoTQo(i-o)’ where aco is the CO-to-molecular gas mass 
conversion factor including He mass. We adopted a conver¬ 
sion factor of aco = 0.8 M© (K km s“^ pc^)“^, the stan¬ 
dard value for ULIRGs (Downes & Solomon 1998). The 
FIR luminosity and dust mass were derived from Lfir = 


dvrAfcJust J'q (z^rest)-S(^rest 5 Thust )d^', aud frOm -i\J^dust — 
/ [( 1 T 2^)z^'cz(z^rest)-S(^rest 5 Thust)] (De Breuck et al. 
2003), where Kd(i^rest) is the dust mass absorption coefficient, 
z/rest is the rest-frame frequency, Tdust is the dust temper¬ 
ature, H(z/rest, Tdust) is the Planck blackbody function, and 
Sohs is the observed flux density. We assume that the ab¬ 
sorption coefficient varies as Kd oc and the emissivity in¬ 
dex lies between 1 and 2 (e.g., Hildebrand 1983), and adopt 
zccz( 125 fim) — 2.64 zb 0.29 m^ kg“^ (Dunne et al. 2003) and 
/3 = 1.5. The dust temperature of Tdust = 39.3 K was derived 
by Negrello et al. (2014) for the integrated value of SDP81, 
which is comparable to high-redshift ULIRGs and SMGs (e.g., 
Kovacs et al. 2006; Symeonidis et al. 2013). We adopted the 
dust temperature for each clump. The SFR was derived from 
SFR = 1.72 X 10 “^°Tfir (Kennicutt 1998a). The derived sur¬ 
face densities of molecular gas and SFR are shown in Table 1. 

5 Discussion 

5.1 Star-formation Properties 

The surface density of molecular gas and SFR for the molec¬ 
ular clumps are plotted in Figure 2. The magnification varies 
with each clump, and a clump with a larger magnification 
traces an intrinsically fainter one. We reiterate that the mag¬ 
nification correction is not necessary when discussing the sur¬ 
face density of molecular gas and SFR, because the gravita¬ 
tional lensing conserves the surface brightness for the identi¬ 
fied clumps. The molecular clumps of SDP81 traced in this 
study have surface densities that are more than three orders 
of magnitude higher than those of local spiral galaxies, and 
they are similar to those of the ~100 pc scale clumps of the 
spatially-resolved SMG SMM J21352 (Thomson et al. 2015). 
We note that the four clumps with high magnification (>100; 
Section 4.1) which are not presented in Figure 2 share the same 
region with the other clumps. The SFR density ranges between 
100-300 M© yr“^ kpc“^, which is lower than the Eddington- 
limited star formation (maximum starburst; Elmegreen 1999) 
(Rybak et al. 2015). It is known that local star-forming galaxies 
follow a relation in the Srs-^^sfr plane (Schmidt-Kennicutt 
relation; Schmidt 1959; Kennicutt 1998a), and it is extended to 
include higher redshift sources (e.g., Genzel et al. 2010). The 
clumps of SDP81 are located in the ‘burst’ sequence (Daddi et 
al. 2010) with a gas depletion time scale of ~ 100 Myr, suggest¬ 
ing active star-forming activity in the clumps. If we adopt the 
Galactic conversion factor (aco = 4.3M© (K km s“^ pc^)“^; 
Bolatto et al. 2013), these clumps are located on the sequence 
of normal star-forming galaxies. Note that if the filling factor 
is less than unity, the data points in Figure 2 will shift to the 
upper-right direction. Regardless of the exact location of these 
clumps in the Ers-^^sfr plane, the derived Esfr/Ers ratios 
are consistent with the values found in the Milky Way clumps 
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Table 1. Properties of molecular clumps. 


ID 

RAirnage 

DeCimage 

RA ^ 

^■^source 

Ds^source 

A^mean 


A^max 

"^mean 

(km s“^) 

Shj" 

(lO'‘M0pc-=) 

^SFR-^ 

(10^ Mq yr “^ kpc “^) 

1 

09:03:11.628 

+00:39:7.63 

09:03:11.558 

+00:39:06.47 

4 

3.6 

4.7 

-94 

1.11 + 0.33 

3.80 + 0.19 

2 

09:03:11.631 

+00:39:7.56 

09:03:11.556 

+00:39:06.47 

4 

3.7 

4.7 

-165 

1.19 + 0.36 

1.60 + 0.20 

3 

09:03:11.636 

+00:39:7.62 

09:03:11.559 

+00:39:06.53 

6 

4.9 

6.6 

-78 

1.13 + 0.35 

1.63 + 0.20 

4 

09:03:11.662 

+00:39:7.25 

09:03:11.560 

+00:39:06.62 

>100 

17 

>100 

118 

1.09 + 0.27 

1.72 + 0.15 

5 

09:03:11.671 

+00:39:7.02 

09:03:11.561 

+00:39:06.62 

>100 

14 

>100 

107 

1.04 + 0.19 

1.58 + 0.11 

6 

09:03:11.672 

+00:39:6.78 

09:03:11.559 

+00:39:06.60 

26 

15 

67 

92 

1.12 + 0.28 

1.79 + 0.16 

7 

09:03:11.661 

+00:39:5.79 

09:03:11.560 

+00:39:06.53 

>100 

17 

>100 

20 

1.10 + 0.36 

1.53 + 0.20 

8 

09:03:11.653 

+00:39:5.57 

09:03:11.562 

+00:39:06.50 

>100 

12 

>100 

-21 

1.05 + 0.16 

1.19 + 0.09 

9 

09:03:11.635 

+00:39:5.39 

09:03:11.560 

+00:39:06.54 

14 

9.7 

21 

40 

1.12 + 0.25 

1.39 + 0.14 

10 

09:03:11.627 

+00:39:5.39 

09:03:11.558 

+00:39:06.60 

8 

6.5 

11 

29 

1.08 + 0.31 

2.91 + 0.17 

11 

09:03:11.444 

+00:39:6.77 

09:03:11.558 

+00:39:06.57 

4 

3.9 

4.4 

92 

1.04 + 0.35 

3.09 + 0.20 

12 

09:03:11.442 

+00:39:6.44 

09:03:11.557 

+00:39:06.49 

4 

3.8 

4.3 

0.9 

1.15 + 0.27 

2.94 + 0.15 

13 

09:03:11.445 

+00:39:6.36 

09:03:11.559 

+00:39:06.47 

4 

4.0 

4.4 

-90 

1.16 + 0.31 

2.08 + 0.18 

14 

09:03:11.448 

+00:39:5.99 

09:03:11.558 

+00:39:06.37 

4 

4.1 

4.6 

-98 

1.13 + 0.32 

1.39 + 0.18 


Note: ID 1 and 13, and ID 4 and 5 are possibly counter images of the same source. “ Right ascension and declination of the ALMA peak positions in the image plane 
( J 2000); ^ Right ascension and declination in the source plane ( J 2000); Mean, minimum, and maximum magnification factors measured within the region of each clump; ^ 
Mean velocity measured within the region of each clump; ® Molecular gas surface density; ^ SFR surface density. 



Ih 2 (Mo pc-2) 

Fig. 2. Molecular gas surface density-SFR surface density plot. The molec¬ 
ular clumps of SDP.81 are plotted as red circles. The vertical and horizontal 
dotted line represents the 2a detection limit on the surface density of molec¬ 
ular gas and SFR, respectively, for a clump with a source size comparable 
to the beamsize. For comparison we also plot data points taken from lit¬ 
erature: Milky Way clouds (Fleiderman et al. 2010; Evans et al. 2014), gi¬ 
ant Hii regions in M33 (Miura et al. 2014), local disk galaxies (Bigiel et al. 
2010; Kennicutt 1998b), local starbursts (Kennicutt 1998b), 2 ; ~ 1-2 star¬ 
forming galaxies (Tacconi et al. 2008), SMGs (Bothwell et al. 2010), and a 
spatially-resolved lensed SMG SMM J21352 (Thomson et al. 2015). The 
solid, dashed, and dot-dashed lines represent the gas depletion time scale 
of 10 Myr, 100 Myr, and 1 Gyr, respectively. If the filling factor is less than 
unity, the data points of SDR81 shift to the upper-right direction. 

(e.g., Heiderman et al. 2010; Evans et al. 2014), suggesting that 
the star-formation law can be applied to molecular clumps at 

z 3. 

5.2 The Nature of SDR81 

By using our lens model (Tamura et al. 2015) and the GLAFIC 
software, we determined the positions of molecular clumps on 


the source plane (Figure 3). The clumps are distributed over a 
2 kpc area, and they overlap with the dust clumps identified in 
Tamura et al. (2015). We also derived the source-plane posi¬ 
tions for the brighter points of the stellar components shown as 
crosses in Figure Id. The points are selected so that the relative 
correspondence of the stellar components between the image 
plane and the source plane is easily visualized. The stellar com¬ 
ponents traced in this study are distributed over kpc in the 
NE-SW direction, which is consistent with the result of Dye et 
al. (2014), and offset from the molecular gas and dust clumps 
toward the west. Since the lensing model of Tamura et al. (2015) 
includes the contribution from the stellar emission detected in 
the HST image, we argue that the extension of the stellar emis¬ 
sion is real and not due to the lens model being based on the 
clumps of AEMA images. We note that, since the main aim of 
this study is to check the positional correlation between molec¬ 
ular gas/dust and stars, we only computed the position of the 
components marked in Figure Id, and the spatial extent of each 
component is not taken into account. 

The mean velocity of each molecular clump measured in the 
image plane is shown in color in Figure 3 (top). A velocity 
gradient with a peak-to-peak value of ~300 km s“^ in the north- 
south direction can be seen, which could indicate a rotating gas 
disk. In Figure 3 (bottom), the molecular gas surface density is 
shown in color scale. The molecular gas surface density appears 
to be higher in the central part of the molecular/dust ‘diskf and 
decreasing toward the edge, although the uncertainty is large. 

What is the nature of SDP.81? One possible scenario is that 
the components of molecular gas, dust, and stars are distributed 
in a several-kpc scale rotating disk, and the stellar emission 
is heavily obscured by dust in the central star-forming region. 
This is qualitatively consistent with the z = 4.05 SMG GN20, 
where an offset of rest-frame UV-emitting region by 4 kpc from 
a region with CO and dust emission has been found (lono et 
al. 2006; Hodge et al. 2015). Alternatively, SDP.81 may be a 
merging system, where dusty starbursts occur in the interface 
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Fig. 3. Positions in the source plane for the molecular clumps, the dust 
clumps, and the stellar components marked in Figure 1. Velocity (top) 
and surface density (bottom) of the molecular clumps are indicated by color 
scale. 


between the two colliding galaxies while the tidal features are 
widely extended. A similar large scale velocity gradient is seen 
in a recent ALMA CO(3-2) image of the Antennae Galaxies 
(Whitmore et al. 2014). The CO(5-4) line we used in this study 
has a relatively high critical density (10^ cm“^) and it may not 
trace the gas of the entire disk. Future observations of the lower 
transition CO lines will allow us to study the detailed distribu¬ 
tion and kinematics of the diffuse gas, providing further insights 
to the exact nature of SDP.81. 
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Data analysis were in part carried out on common use data analysis 
computer system at the Astronomy Data Center, ADC, of the National 
Astronomical Observatory of Japan. 
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